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Abstract 

We have studied spin-spin correlation functions in the ordered 
phase of the two-dimensional g-state Potts model with q = 10, 15, 
and 20 at the first-order transition point (3 t - Through extensive Monte 
Carlo simulations we obtain strong numerical evidence that the cor- 
relation length in the ordered phase agrees with the exactly known 
and recently numerically confirmed correlation length in the disordered 
phase: (, {Pt) = £,d{ftt)- As a byproduct we find the energy moments 
in the ordered phase at fit in very good agreement with a recent large 
(/-expansion. 



PACS numbers: 05.50.+q, 75.10.Hk, 64.60.Cn, 11.15.Ha 



1 Introduction 



First-order phase transitions have been the subject of increasing interest in 
recent years. They play an important role in many fields of physics as is 
witnessed by such diverse phenomena as ordinary melting, the quark decon- 
finement transition or various stages in the evolution of the early universe [1]. 
Even though there exists already a vast literature on this subject [2], many 
properties of first-order phase transitions still remain to be investigated in 
detail. Examples are finite-size scaling (FSS) [3], the shape of energy or 
magnetization distributions [4, 5], partition function zeros [6], etc., which 
are all closely interrelated. An important approach to attack these problems 
are computer simulations [7]. Here the available system sizes are necessarily 
limited and reliable FSS analyses are of utmost importance. The distinguish- 
ing input parameter in such analyses is the correlation length £ which sets 
the relevant length scale of the system and thus the range of validity of the 
commonly employed FSS Ansatze [3]. 

The well-known paradigm to investigate these questions in detail is the q- 
state Potts model [8] which can be tuned from a second-order through weakly 
first-order to a strong first-order transition by varying the number of states 
q. In two dimensions many quantities of interest are known exactly [9] and 
have been shown to have interesting algebraic interpretations in knot theory 
[10]. The theoretical knowledge of correlation lengths, however, is still quite 
limited. Only in the disordered phase exactly at the first-order transition 
point j3 t (q > 5) an explicit formula for the correlation length £d(A) is known 
[11, 12]. By analogy with the Ising model [13] and on the basis of previous 
numerical data [14] it has been speculated [12] that the correlation length 
in the ordered phase is given by £ (A) = ^d(Pt)- In f ac t, very recently this 
could be proven for a certain definition £ 0jl , but the possibility of a larger cor- 
relation length £ 0j2 > £ ,i was left open [15]. In this note we report a careful 
Monte Carlo study of the leading correlation length in the two phases which 
yields unambiguous numerical evidence against this conjecture. Rather, our 
numerical analysis strongly favors an earlier conjecture that £ (/3t) = 
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2 Model and simulation 



The Potts model is defined by the partition function 

z = E e-/3£ ; E = - E<W Si = i,..., q , (i) 

where i denote the lattice sites, (ij) are nearest-neighbor pairs, and S SiSj is 
the Kronecker delta symbol. We simulated the model with q = 10, 15, and 20 
in the ordered phase at j3 t = ln(l + s/q) on 2D lattices of size V = L x x L y . 
For each q we considered two different lattice geometries, namely L x L 
and 2L x L lattices with L = 150, 60, and 40 for q = 10, 15, and 20, 
respectively. To take advantage of translational invariance we used periodic 
boundary conditions. From our experience with simulations in the disordered 
phase we knew that lattice sizes L m 14£ d are large enough to suppress 
tunneling events. In fact, starting from a completely ordered configuration, 
we never observed a tunneling event into the disordered phase, thus allowing 
statistically meaningful pure-phase measurements of energy moments and 
correlation functions. In Fig. 1 we show the probability distributions P(e) of 
the energy (normalized to unit area) in the ordered as well as in the disordered 
phase, demonstrating that the two peaks are indeed very well separated. 

While in our earlier study [16] of the disordered phase we employed the 
single-cluster algorithm, here we found it more efficient, in terms of real- 
time performance on a Cray-YMP, to update the spins with a vectorized 
standard heat-bath algorithm. With this update algorithm the value of the 
spins in the largest (spanning) cluster never changed, as would clearly be the 
case with fixed boundary conditions. The main difference between the two 
types of boundary conditions is that, with fixed boundary spins, excitations 
of disordered bubbles would be repelled from the walls, while with periodic 
boundary conditions they can freely move around the torus. The statistic 
parameters of the runs are collected in Table 1, where we give the number of 
update sweeps in units of the integrated autocorrelation time of the energy, 
Ti nt e , which is roughly the same for the two lattice geometries. 

To determine the correlation length we followed Gupta and Irback [14] 
and considered the k y n ^ = 2im/L y momentum projections (i = (i x ,i y )), 

g {n) (i*,j x ) = y E - -)e ik ^^\ (2) 
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For n ^ the projections are free of constant background terms and, similar 
to our analysis in the disordered phase, we tried to determine £ from fits of 

g {n \x) = g (n) (i x ,0) to the Ansatz 

9 (n \x) = ach(^|^) + 6ch(c^|^), (3) 

with f W = £ /^l + {27rn£ /L y ) 2 . Tests for the Ising model suggested that 
this is a good approximation for L/n£ > 10. For the Potts models with 
L/£ ~ 14 we expect that ffl is smaller than by about 10% for the k y - 
projections and 2% for the /^-projection on the 2L x L lattice, respectively. 
Even though in the ordered phase it is more efficient to update the spins 
with a heat-bath algorithm, for the measurements we decomposed the spin 
configuration into stochastic Swendsen-Wang cluster and used the cluster 
estimator (S SiSj - l/q) = (1 - l/q)(Q(i,j)), where Q(i,j) = 1, if i and j 
belong to the same cluster, and 6 = otherwise. All error bars are estimated 
by means of the jack-knife technique. 



3 Results 

We first checked that the average energy agrees with the exact result [9] , and 
compared the second and third energy moments, c Q = = (3 2 V{{e — 

(e)) 2 ) and ^ = V 2 ((e — (e)) 3 ), with the large q expansions of Ref.[5]; cp. 
Table 2. Even for our smallest value q = 10 and the third moment we observe 
a very good agreement with the (Pade resummed) large q expansion. In addi- 
tion we looked at the magnetization m = (q (max{nj})/V — l)/(q— 1) and its 
cluster estimator m! = (|C| span )/V, where denotes the number of spins of 
"orientation" i = 1, . . . , q and |C| span is the size of the largest (spanning) clus- 
ter. Table 3 shows that the two estimators give almost identical results which 
agree very well with the exact expression m = Yi^Li [(1 — x n )/(l — x 4n )], 
where, for q > 5, x is defined by q = x + 2 + x~ x (0 < x < 1) [17]. Also 
shown is the susceptibility 1 , Xo = = V{(m — (m)) 2 ), and the third 

moment, mff 1 = V 2 ((m — (m)) 3 ). 

x To conform with the present normalization, our Monte Carlo estimate for \d [16] 
should be multiplied by a factor of l/(q — 1). 



3 



For the estimate of we concentrated on the k y = 27r/L y projec- 

tion g^(x). The qualitative behaviour of is illustrated in the semi-log 
plots of Fig. 2. For comparison also our previous results [16] for in 
the disordered phase are shown. Already these plots suggest that the two 
correlation functions are governed by the same asymptotic decay law, i.e., 
that £ (Pt) = €d(Pt)- In f ac t, the dotted lines interpolating the g^ data are 
constrained fits to the Ansatz (3) assuming that £ = ^ (= 10.559519..., 
4.180954 . . ., 2.695502 ... for g = 10, 15, 20 [11, 12]). To be sure we also per- 
formed unconstrained four-parameter fits using Ansatz (3). Figure 2 shows 
that, despite our high statistics, in the ordered phase it was impossible to 
get reliable estimates of g^ for very large distances. Compared with our 
analysis in the disordered phase the available fit intervals are consequently 
shifted to smaller x. This makes the fits in the ordered phase somewhat more 
sensitive to higher-order excitations and, based on our previous experience 
with varying fit intervals, we estimate that the numbers for ffl collected in 
Table 4 should underestimate the true value by about 15 — 30%. For a com- 
parison of with £d(Pt) we have therefore used approximately the same 
fit intervals. For example, for q = 10, L x L lattice and the fit starting at 
^min = 11, we obtain from Table 4 £ G (/5 t ) = 7.9(9). Recalling our estimate of 
id(fit) = 8.8(3) for a comparable fit in the disordered phase [16], this yields 
a ratio of £ (Pt) /€d(Pt) = 0.9(1), suggesting again that = This 
is corroborated by the corresponding estimates for q = 15 and q = 20, which 
are actually closer to unity. 

To make this statement even more convincing we have plotted in Fig. 3 the 
ratio £® ff /Q s , where £ eff = 1/ \n[g(x)/g(x+l)] is the usual effective correlation 
length (with the correction £ = £ (£o^) already taken into account). Here 
we have again implicitly assumed that higher-order excitations play a similar 
role in both phases and effectively drop out when plotting the ratio of the two 
correlation lengths. Figure 3 shows that, over a sizeable range of distances, 
£o ff /£d ff = 1) with an uncertainty of only about 5%. We interpret this as 
further support of our conjecture that £ (/5t) = £d(/3t)- 

4 Discussion 

To summarize, applying a non-zero momentum projection to the correlation 
function in the ordered phase we determined £ {f3t) with an accuracy compa- 
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rable to previous investigations in the disordered phase [16]. By comparing 
the two correlation lengths at fl t we obtain strong numerical evidence that 
£o = £d- At first sight this is in striking disagreement with a very recent 
exact proof [15] of the earlier conjecture £ = for one definition of the or- 
dered correlation length, £ o l . For another definition, £ o 2 , however, only the 
relation £ o 2 > £ 0j i could be established in Ref.[15]. This is clearly consistent 
with our result if we identify the numerically determined £ with £ 0j 2. We 
are currently investigating this problem in more detail [18] by using precisely 
the definitions of Ref.[15] which are based on geometrical properties of Potts 
model clusters such as, e.g., the distribution function of the cluster diameter. 
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Figure 1: Probability distribution P(e) in the two phases for L x L (dotted 
line) and 2L x L lattices. 

Table 1: Integrated autocorrelation time Ti ntje of the energy and the number 
of update sweeps in units of Ti ntj e- 





q = 10 


q = 15 


q = 20 




w 140 


w 20 


w 10 


Lx L 
2L x L 


118 000 
121000 


1280 000 
1280 000 


1280 000 
5120 000 



Table 2: Comparison of numerical and analytical results for energy moments 
at j3 t in the ordered phase. 



Observable 


g = 10 


q = 15 


q = 20 


e (MC, Lx L) 
e (MC, 2L x L) 
e Q (exact) 


-1.664177(81) 
-1.664262(57) 
-1.664253... 


-1.765850(34) 
-1.765875(26) 
-1.765906... 


-1.820722(43) 
-1.820689(14) 
-1.820684... 


c (MC, Lx L) 
c (MC, 2L x L) 
c (large q) 


17.95(13) 
17.81(10) 
18.1(1) 


8.016(21) 
8.004(19) 
8.00(3) 


5.351(15) 

5.3612(55) 

5.362(3) 


A^(MC, Lx L) 
//( 3 ) (MC, 2L x L) 
//^ 3 ) (large g) 


1979(87) 
1836(71) 
1850(40) 


180.5(3.1) 
189.7(5.1) 

179(4) 


57.0(1.3) 
56.24(40) 
56.8(4) 
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Table 3: Comparison of numerical and analytical results for magnetization 
moments at (3 t in the ordered phase. 



Observable 


q = 10 


q = 15 


q = 20 


m (MC, Lx L) 
m' (MC, Lx L) 
m (MC, 2L x L) 
m' (MC, 2L x L) 
m (exact) 


0.857047(71) 
0.857047(71) 
0.857113(49) 
0.857113(49) 
0.857106... 


0.916631(21) 
0.916634(21) 
0.916648(16) 
0.916648(16) 
0.916663... 


0.941199(21) 

0.941197(21) 

0.9411782(66) 

0.9411791(66) 

0.9411759... 


Xo (MC, Lx L) 
Xo (MC, 2L x L) 


4.750(60) 
4.663(43) 


0.8090(36) 
0.8095(38) 


0.3348(17) 
0.33509(55) 


(MC, Lx L) 
mf (MC, 2L x L) 


-1521(85) 
-1372(62) 


-45.9(1.2) 
-49.4(2.2) 


-8.55(32) 
-8.321(88) 



Table 4: Numerical estimates of the correlation length £^(A) from four- 
parameter fits to the Ansatz (3) in the range x m i Q — x max . For the 2L x L 
lattices the fits along the x and y direction are distinguished by the index. 
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= 40 


Lattice 


x . - 

mm 


-x 

max 




x . - 

mm 


-x 

max 




x . - 

mm 


-x 

max 






4- 


-54 


7.4(4) 


3- 


19 


3.0(1) 


2- 


12 


2.0(1) 


Lx L 


7- 


-54 


7.5(6) 


5- 


19 


3.2(3) 


3- 


12 


2.0(2) 




11- 


-54 


7.5(8) 


7- 


19 


3.4(9) 


4- 


12 


2.0(2) 




4- 


-40 
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3- 


22 


3.2(1) 
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14 
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7- 


-40 
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5- 


22 
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14 
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2.1(1) 


(2L x L) y 


7- 


-40 
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22 
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3- 


14 


2.2(1) 
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22 
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14 


2.2(2) 
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Figure 2: Semi-logarithmic plots of the projected correlation function 
for q = 10, 15 and 20 in the ordered phase. For comparison also in the 
disordered phase [16] is shown. 
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Figure 3: Ratio of effective correlation lengths in the ordered and disordered 
phase for q = 10, 15 and 20. 
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